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ABSTRACT 

It has been argued that the flux anomalies detected in gravitationally lensed QSOs are 
evidence for substructures in the foreground lensing haloes. In this paper we investigate 
this issue in greater detail focusing on the cusp relation which corresponds to images of 
a source located to the cusp of the inner caustic curve. We use numerical simulations 
combined with a Monte Carlo approach to study the effects of the expected power law 
distribution of substructures within ACDM haloes on the multiple images. 

Generally, the high number of anomalous flux ratios in the cusp configurations is 
unlikely explained by 'simple' perturbers (subhaloes) inside the lensing galaxy, either 
modeled by point masses or extended NFW subhaloes. We considered in our analysis 
a mass range of 10 5 — 1O 7 M0 for the subhaloes. We also demonstrate that including 
the effects of the surrounding mass distribution, such as other galaxies close to the 
primary lens, does not change the results. We conclude that triple images of lensed 
QSOs do not show any direct evidence for dark dwarf galaxies such as cold dark matter 
substructure. 
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1 INTRODUCTION 



Cold Dark Matter (CDM) simulations predict many more 
low mass satellite haloes than are actually observed in the 
Milky Way (Klypin et al. 1999, Moore et al. 1999). It seems 
■ that 10-15% of the mass was left in satellites with perhaps 
1-2% at the projected separations of 1-2 Einstein radii (Re) 
where we see most lensed images (e.g. Zentner & Bullock 
2003, Mao et al. 2004); this is far larger than the observed 
fraction of 0.01-0.1% in observed satellites (e.g. Chiba 2002). 
Solutions to this mismatch were proposed in three broad 
classes: satellites are present but dark if star formation is 
prevented (Bullock, Weinberg & Kravtsov 2000), satellites 
are destroyed due to self-interacting dark matter (Spergel & 
Steinhardt 2000), or their formation is prevented by chang- 
ing the power spectrum to something similar to warm dark 
matter with significantly less power on the relevant mass 
scales (e.g. Bode et al. 2001). These hypotheses left the ma- 
jor observational challenge of distinguishing dark satellites 
from non-existent ones. This became known as the CDM 
substructure problem. 

It has been argued that a possible signature of the pres- 
ence of dark matter substructures can be found in strong 
gravitational lensing of QSOs (Mao & Schneider 1998; Met- 
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calf & Madau 2001; Chiba 2002; Metcalf & Zhao 2002; Dalai 
& Kochanek 2002: Kochanek & Dalai 2004). If a distant im- 
age source is close to a cusp (from inside) in a caustic curve, 
three of the images will be clustered together and the sum 
of their magnifications will be zero (Zakharov (1995), taking 
the negative parity image to have negative magnification). 
This relation holds for a wide class of smooth analytic lens 
models (Keeton et al. 2003); on the other hand all known 
observed lensed QSOs violate this relation. This has been 
explained with the presence of cold dark matter substruc- 
tures within the lensing galaxy's halo. 

However, the discrepancy found in some systems may 
be due to microlensed stars rather than to cold dark matter 
substructures (Keeton et al. 2003), even if the most peculiar 
problem is the anomalous flux ratios in radio lenses. Ra- 
dio sources are essentially unaffected by the ISM of the lens 
galaxy (see however Koopmans et al. 2003), true absorption 
appears to be rare, radio sources generally show little vari- 
ability and most of the flux should come from regions too 
large to be affected by microlensing. Therefore dark matter 
subhaloes appear to be the most likely explanation. 

By using low resolution simulation of galaxy formation 
Bradac et al. (2004) claimed that the level of substructures 
present in simulation produces violations of the cusp rela- 
tion comparable to those observed. Amara et al. (2004) im- 
planted an idealized model of a galaxy into the center of a 
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high resolution galactic halo extracted from dissipationless 
N-Body simulations to test the effects of substructures on 
lensed images. Their findings contrast those of Bradac et al. 
(2004), since they found that the substructures produced in 
a ACDM halo are not abundant enough to account for the 
observed cusp caustic violation. The results of Amara et al. 
(2004) were also confirmed in a recent work by Maccio et 
al. (2005). In the later work, in which a fully hydrodynami- 
cal simulation of galaxy formation is used, it is shown that 
the presence of a dissipative component greatly enhances 
the surviving probability of satellites, expecially close to the 
center of the galaxy. Nevertheless Maccio et al. (2005) also 
demonstrated that the impact on lensing of subhaloes in the 
mass range 10 7 — 10 10 M© is very small. Even with a number 
of subhaloes about 8 higher of the observed one in this mass 
range, the number of multiple lensed QSOs that show a vio- 
lation of the cups relation is less than 24%, in contrast with 
an observed one of about 60%. This means that if the vio- 
lation of the cusp relation is due to substructures inside the 
primary lens, these must have a mass smaller than 10 7 Mq. 

The aim of this work is to study the influence of sub- 
haloes with mass 10 5 — 1O 7 M0 on the cusp relation violation. 

The outline of the paper is the following: first we briefly 
summarize the cusp relation, then in section 3 we present the 
lensing numerical simulations and our modeling for the pri- 
mary lens, subhaloes and extra-haloes. Tests for our models 
and results for three different cusp configurations are pre- 
sented in section 4. Section 5 is devoted to a short discussion 
on the fold relation. A discussion of the results and our con- 
clusions are presented in section 6. Throughout this paper 
the single large halo that is causing the QSO multiple images 
is referred to as the primary lens. The additional small scale 
haloes (inside the host halo) are referred to as subhaloes or 
substructures. Haloes beyond the virial radius of the pri- 
mary lens are referred to as extragalactic haloes. We adopt 
the standard ACDM cosmological model with the following 
parameters Q m = 0.3, = 0.7, as = 0.9 and H — 70 km 
s" 1 Mpc- 1 . 



2 THE CUSP RELATION 

There are basically three configurations of four-image sys- 
tems: fold, cusp, and cross (Schneider & Weiss 1992). In this 
paper we will mainly concentrate on the cusp configuration, 
that corresponds to a source located close to the cusp of the 
inner caustic curve. The behavior of gravitational lens map- 
ping near a cusp was first studied by Blandford & Narayan 
(1986), Schneider & Weiss (1992), Mao (1992) and Zakharov 
(1995), who investigated the magnification properties of the 
cusp images and concluded that the sum of the signed mag- 
nification factors of the three merging images approaches 
zero as the source moves towards the cusp. In other words 
(e.g. Zakharov 1995) : 



\^a\ + |ms| + |mc| 



0, for fitot — » oo (1) 



where p tot is the unsigned sum of magnifications of all four 
images, and A,B & C are the triplet of images forming the 
smallest opening angle (see figure 7). By opening angle, we 
mean the angle measured from the galaxy center and being 
spanned by two images of equal parity. The third image lies 



inside such an angle. This relation is an asymptotic relation 
and holds when the source approaches the cusp from inside 
the inner caustic "astroid" . This can be shown by expanding 
the lensing map to third order in the angular separation from 
a cusp (Schneider & Weiss 1992). Small scale structure on 
scales smaller than the image separation will cause i? CU sp to 
differ from zero fairly independently of the form of the rest 
of the lens. Indeed, a substructure is more likely to reduce 
the absolute magnification for negative magnification images 
(Metcalf & Madau 2001, Schechter & Wambsganss 2002, 
Keeton et al. 2003) and to increase it for positive parity 
images. 



3 LENSING SIMULATIONS 

We use the lensmodel package (Keeton 2001)^ modeling the 
main lens galaxy as a singular isothermal ellipsoid (SIE) 
and the substructures as NFW (Navarro, Frenk & White 
1997) haloes. First, using the gravlens task, we find three 
lens configurations for which the cusp relation is roughly 
satisfied (figures 7, 10 and 16). 

As second step a variable number of substructures is 
added to the main lens (see section 3.2 for details on their 
number density and physical properties). For this new lens- 
ing system (main lens plus subhaloes) we compute again 
positions and fluxes of the images (subhaloes mainly tend 
to modify fluxes more than positions, see Kochanek (2004) 
and section 4.1), obtaining a new value for the cusp relation 
i?cus P . This procedure is repeated more than 20.000 times 
for each of three studied positions of the source (figures 7, 
10 and 16): this allows us to compute the probability distri- 
bution of the .Rcusp value in presence of subhaloes (i.e. figure 



3.1 Primary Lens 

The observed discrepancy in the flux ratios, compared with 
the expected universal relation from a cusp or fold singular- 
ity, suggests that it is an intrinsic difficulty for smooth lens 
models, not associated with a particular parameterization. 
For the scope of this paper it is sufficient to choose just a 
single smooth lens model for the primary lens. Therefore, 
we select, as a smooth lens model, a singular isothermal el- 
lipsoid (SIE) (Kormann, Schneider, & Bartelmann 1994) to 
take advantage of its simplicity. This model has been widely 
used in lens modeling and successfully reproduces many lens 
systems (e.g. Keeton et al. 1998, Chiba 2002, Treu & Koop- 
mans 2004). An isothermal profile for the total mass distri- 
bution of elliptical galaxies is well supported by the detailed 
dynamical studies of local ellipticals (Gerhard et al. 2001), 
individual lens modeling, and statistics (e.g. Maoz & Rix 
1993; Kochanek 1995; Grogin & Narayan 1996). 

The ellipsoidal primary lens has a mass equal to 5 x 
1O 11 M0, it is oriented with the major axis along the y axis 
in the lens plane and has an ellipticity of 0.33. The redshifts 
of the lens and he source are fixed to zi = 0.3 and z s — 1.71 



t The software is public available via web site: http://cfa- 
www.harvard.edu/castles 
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respectively, agreement with the typical observed ones (in 
this case we use PG1 115+080 data, see Tonry 1998) 



3.2 Subhaloes 

Since it has been shown that the number density of sub- 
haloes with mass M > 10 7 Mq is not sufficient to explain the 
observed number of violation in the cusp relation (Amara et 
al. 2004, Maccio et al. 2005) , the aim of this work is to inves- 
tigate the impact of substructures below this mass threshold 
that is fixed by the current resolution limits of numerical 
simulations. 

We would like to emphasize that we are using only one 
lensing plane, this means that we will consider only effects 
due to substructures being at the same redshift of the main 
lens (see Chen, Kravtsov & Keeton 2003, Metcalf 2005 for 
an estimation of the effects of haloes along the line of sight). 
In order to evaluate the number density substructures in the 
mass range 10 5 - 10 7 M Q we have made some extrapolations 
based on results from high resolution N-body simulations. 
The mass function of subhaloes inside the virial radius of an 
halo is close to a power law (Diemand et al. 2004, DMS04 
hereafter, Gao et al. 2004, Reed et al. 2005): 

N(> m) oc m~ 9 , (2) 

with a slope (3 w 1, so that we expect to have a factor » 100 
more subhaloes inside the viral radius if we move our mass 
threshold from 10 7 M Q to 1O 5 M . 

As said in the previous section such small haloes will af- 
fect the -Rcusp relation only if their distance from the images 
is of the same order or smaller than the distance between 
the images themselves. Therefore we need an estimation of 
the number of haloes inside a small area surrounding the 
images. This number will also depend on the distance of our 
area from the center, due to the fact that the number den- 
sity of haloes increases approaching the center of the main 
halo (primary lens) as clearly shown in fig 1, which is based 
on numerical simulations of 4 galaxy size haloes (DMS04) . 

Consequently the number of subhaloes with a mass 
greater than m inside an area A at a distance R from the 
center of the galaxy is: 

(N rv (>m ))^N(R)A 



N A (>m,R) 



(3) 



where (N rv (> mo)) is the average number density of sub- 
haloes with m > mo (being mo an arbitrary mass value) 
inside the virial radius r v and N(R) is the radial 2D num- 
ber density of satellites at a projected distance R from the 
center in units of {N rv (> mo)) (see figure 1). These last 
two quantities can be obtained directly from N-body simula- 
tions. In the following we will use results from DMS04 (table 
1 of their paper, simulations G0-G3). Maccio et al. (2005) 
have shown that the presence of baryons inside subhaloes 
enhance the probability to find haloes close to the center of 
the galaxy with respect to results from dissipationless simu- 
lations. This is true for satellites with m > 5 x 10 7 Mq , which 
are massive enough to retain baryons inside their potential 
well and then form stars. Since we do not expect such effect 
for the mass scales involved in this work (« 10 5 - 1O 6 M ) 
we can use pure N-body simulation as starting point for our 
analysis. 

Since the typical separation between images in lensed 



10 
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Figure 1. Two dimensional radial number density of subhaloes 
in units of the average number density inside the virial radius. 
Result obtained averaging the 4 high resolution galaxies presented 
in DMS04, using three different projections for each galaxy. 



QSOs is roughly a few arcsec, we fix R « 1 arcsec and A = 6 
arcsec 2 {Aq hereafter). We remind that at a redshift of zi — 
0.3 one arcsec corresponds to 4.55 kpc for the cosmological 
model adopted in this paper. Using eq. 3 and adopting a 
mass threshold for substructures of m — 5 x 1O 5 M0 the 
number of subhaloes inside A ranges from 4 to 12 (depending 
on the uncertainties on N(R), see figure 1). For two 1O 7 M0 
haloes the surface mass density within the selected area A 
is O.69/i _2 M pc~ 2 , this means a fraction w 10~ 3 of the 
total dark matter surface density in substructures, in good 
agreement with the results of Mao et al 2004. 

For each lensing configuration analyzed in this work we 
added a random number of substructures between 4 — 12 to 
the primary lens with a random mass generated according to 
cq. 2 in the range 5 x 10 5 — 1O 7 M0. These subhaloes are then 
placed following the 2D density profile inside the area A& 
that encloses the three images (cfr the (blue) square in figure 
7). We have modeled our subhaloes with an NFW (Navarro, 
Frenk & White 1997) density profile; for the « 10 6 M Q sub- 
haloes relevant for lensing substructure studied in this work, 
the NFW profile inferred from N-body simulation is the 
most natural choice, because on these mass scales the ef- 
fect of baryons (that are able to modify the slope of the 
density profile for greater masses (Maccio et al. 2005)) is 
very tiny because the potential well of these haloes is not 
deep enough to retain them expecially in presence of a ion- 
izing background. We have adopted different concentration 
parameters (see sec. 4.1) to mimic the scatter present in the 
mass-concentration relation (Bullock et al. 2001). 



3.3 Extrahaloes 

Lensing galaxies are not isolated object, since they usually 
belong to group of galaxies (Keeton et al. 2000). Moreover 
each galaxy has its own satellites galaxies with masses in 
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the range 10 9 - 10 n M Q . Consequently we have modeled the 
presence of these extrahaloes in the same way of the inner- 
haloes. In this work we call substructures or innerhaloes, 
haloes with mass < 1O 7 M0 that are close to the image 
position; we reserve the term extrahaloes for haloes with 
M > 10 9 Mq. We do not consider extrahaloes that for pro- 
jection effects can lie inside the primary lens close to the 
images positions (i.e. Oguri 2005). 

We considered three different categories of extrahaloes: 
i) haloes with mass 10 9 < M < 10 10 M Q and with a pro- 
jected distance r between 60 and 200 kpc, these repre- 
sent the satellites galaxy of the primary lens (the expected 
number for these haloes can be estimated again form N- 
body simulation and it is roughly 6-8); ii) haloes with mass 
10 11 < M < 10 12 M Q and distance 300 < r < 700 kpc in or- 
der to mimic the presence of companion galaxies; iii) haloes 
with 10 12 < M < 5 X 10 13 M Q and 700 < r < 1200 kpc to 
take into account the possible presence of a nearby cluster 
of galaxies. 

The number of extrahaloes has been fixed between 2-8 
and 2-4 for the ii) and iii) case respectively, as suggested 
by observations/simulations (Metcalf 2005, Amara 2004). 
While the extrahaloes in case i) are placed in a circularly 
symmetric way around the center of the galaxy, the posi- 
tion of groups and clusters of galaxies must be modeled in 
an asymmetric way. Therefore we placed them only in the 
quadrant with positive coordinates in the lens plane (being 
the lens in [0 : 0]). We used SIE as lens model for extrahaloes 
to take into account the presence of baryons their- inside and 
we have generated 20.000 different configurations. 



4 RESULTS 

In this section we present results of our Monte Carlo sim- 
ulations. We have analyzed 3 cusp configurations: Configl 
(figure 7), Config2 (figure 10) and Config3 (figurel6). They 
mainly differ for the value of R CU s P in the unperturbed case 
that grows from 0.01 for Configl to 0.243 for Config3, due 
to a different position of the source inside the inner caustic 
curve. 

4.1 Testing our model 

Before proceeding further in our analysis we present some 
tests on the parameter adopted in our lensing simulation. 
The presence of substructures acts both on positions and 
fluxes of the images. When a substructure is added the 
gravlens code adjusts the positions and fluxes of the images 
by minimizing the \ 2 between the old (unperturbed) and 
new (perturbed) positions/fluxes. The value of \ 2 depends 
on the error that we assign to the unperturbed positions 
and fluxes (usually these are the observational errors); this 
means that if the error on positions is smaller than the er- 
ror on fluxes the code will change the latter ones more than 
moving the images to obtain a lower value for \ 2 - 

Therefore in principle the value of the perturbed i? C us P 
is influenced by the error assigned to the position of the 
images, fn figure 2 we clearly show that this effect is very 
small even for a big variation of the image position errors. 
This plot illustrates the probability distribution for 7? cuap 
(obtained using 20.000 realizations) when 4 substructures 




0.04 0.08 0.12 0.16 0.2 

Figure 2. Probability distribution of i? CU s P values for 2 substruc- 
tures inside A3 for different position errors. 

are added to the primary lens : the solid line is for an error 
on positions of ej = 10 -2 arcsec (« 5%) the dashed one for 
£2 = 10 arcsec (i.e. the code has complete freedom in mov- 
ing the images). The two distributions of the i? C us P values 
are very similar, and this also confirms findings of other au- 
thors: the influence of substructures on the image positions 
is less strong than the one on fluxes (Kochanek 2004). We 
will adopt 10~ 2 arcsec for the error on positions for all our 
simulations. 

A correct determination of the subhaloes properties (see 
section 3.2) is a key ingredient in computing the -R CU s P value: 
figure 3 shows the influence of the number density and mass 
range of substructures on the R CU s P relation. Both changing 
the area in which subhaloes are distributed (keeping fixed 
their number: left panel) or their mass range (right panel) 
leads to complete different results. 

As second step we have tested if the subhaloes that 
live outside the small area (A) surrounding the images can 
substantially modify the i? CU s P relation. For this purpose 
we have defined three different areas for substructures for 
Config2: A 3 = [-1.5 : 1.5] x [1 : 2], A 6 = [-1.5 : 1.5] x [0 : 2] 
and A12 = [-1.75 : 1.75] x [-1.25 : 2.25]. The last one (A 12 ) 
is big enough to cover all the image positions. Changing the 
size of the area and his position in respect to the center 
of the lens, the number of substructures changes according 
to eq: 3. The number of subhaloes inside the three areas is 
4, 9 and 19 respectively. Figure 4 shows our results. The 
-Rcus P probability distribution is weakly affected by the size 
of the area. This indicates that are the subhaloes close to the 
image positions the ones responsible for the R cusp relation 
modification. 

Figure 5 shows results for Config2 (figure 10) of the 
.Rcusp probability distribution for 3 subhaloes with mass 
5 x 10 5 < m < 1O 7 M inside A 3 (solid line) and for 8 sub- 
haloes with 10 7 < m < 10 9 Mq inside an area of 57 arcsec 2 
([—5:5] x [0:6] with the exclusion of A3). For this second 
population of substructures both the number density and 
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Figure 3. Probability distribution of i? cusp with 4 substructures: 
left panel influence of the substructure mass. Right panel effects 
of the projected number density. 



Figure 5. Effects on the -Rcusp relation of substructures inside 
a small area (A3) surrounding the images (solid line) and in a 
larger area [—5 : 5] X [0 : 6] (57 arcsec) outside A3 (dashed line). 
In the first case the subs are in the mass range [5 X 10 5 : 10 7 Mq] 
in the second one [10 7 : IO 9 M ]. 
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Figure 4. Probability distribution of i? CU sp for substructures dis- 
tributed on different areas. 

the masses are over estimated by a large fraction. Neverthe- 
less its effect on the cusp relation is very small and the -Rcusp 
probability distribution is close to a delta function centered 
on the unperturbed value. 

Modeling our subhaloes as NFW structures we have one 
more free parameter: the concentration. It is well know that 
the concentration correlates with the mass of the halo, even 
if a consistent scatter is present in this relation. Extrapolat- 
ing results from Bullock et al 2001, the mean concentration 
in our mass range is around 33. In order to test the influ- 
ence of the concentration of our subhaloes in modifying the 



.Rcusp relation we repeated our analysis keeping fix mass and 
position of the subhaloes and varying their concentrations. 
Results are shown in figure 6. As expected denser haloes 
have a stronger impact on the -Rcusp relation. We have also 
considered two other mass profiles for the subhaloes: SIS 
(singular isothermal sphere) and point like approximation. 
The first one, less favored by simulations, can be seen as 
an upper limit for the NFW profile, due to the fact that 
its inner slope for the density profile is proportional to r~ 2 . 
For this kind of subhalo model the Rcusp probability is just 
slightly above the ones for NFW haloes with c=35, so we do 
not expect a big change in our results using SIS instead of 
NFW subhaloes. On the other hand the point mass approx- 
imation leads to an underestimate of the effect of subhaloes 
on JJcusp since for a fixed mass a point-like object has a 
smaller Einstein radius than an NFW halo (Keeton 2003). 
In the following we have adopted a concentration parameter 
of 35 for our subhaloes in order to try to maximize their 
effect on the -R CU sp relation. 

4.2 Configuration 1 

In the first configuration analyzed (Configl hereafter) the 
source is close to the right cusp of the inner caustic curve 
(see fig 7) and the unperturbed -Rcusp value is 0.01. The 
critical and caustic curves refer to the unperturbed case but 
they are not different from perturbed ones at the level of 
resolution. 

For this configuration we have generated 20.000 dif- 
ferent lensing systems that include substructures according 
to eq: 3. In figure 8 is shown the probability distribution 
for .Rcusp for different numbers of substructures. The max- 
imum of the probability is obtained for the unperturbed 
value (0.01) and the tail of the distribution extends to 
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Figure 6. i?, C usp distribution probability for different subhaloes Figure 8. Probability distribution of R-variation for a different 

mass profiles: NFW profile with different concentration values, number of substructures inside Aq. 

SIS profile and point-like approximation. 
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Figure 7. Unperturbed lens configuration: Configl (R=0.01) 



Rcusp = 0.12 but for a very low number of configurations 
(less than 1.0%). Figure 9 shows in a logarithmic plot the 
tail of the distribution presented in figure 8: it is possible 
to note that an increase of the total number of substruc- 
tures (from 4 to 6) does not substantially change the value 

of R cusp • 

From figure 8 one see that in the 10% of the configu- 
rations the final value of -Rcusp is even less than the unper- 
turbed value and it is closer to the theoretical expectation 
of Rcusp = (see next section on the second configuration 
for more details). 
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Figure 9. Same of figure 8 with logarithmic scale for the y axis. 
4.3 Configuration 2 

In this configuration (Config2 hereafter) the source is close 
to the upper cusp of the inner caustic curve (see fig 10) with 
an unperturbed value of Rcusp = 0.09. 

Figure 11 shows the analogous of figure 8 for Config2 
and figure 12 shows the tail of the distribution for large 
values of i? C us P . Even in this case the maximum of the prob- 
ability distribution is centered on the unperturbed value. 
Here the distribution of the values of R CU sp is more symmet- 
ric than for Configl and it is more evident that the effect of 
substructures not only increase the value of the cusp ratio 
but can also reduce it. 
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Figure 10. Second unperturbed lens configuration: Config2 
(R=0.09) 



Figure 12. Same of figure 11 with logarithmic scale for the y 
axis. 
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Figure 11. Probability distribution of R-variation for config 2, 
for a different number of substructures inside A§. 



To better illustrate this effect we have isolated one of the 
configurations in which we find a reduction of R C us P (figure 
13). When the distribution of subhaloes is non symmetric in 
respect to the triplets of images and one of the perturbers is 
close to one of the external images the latter image results 
to be more magnified than the others. In the unperturbed 
configuration |/^(-B)| > fi(A) + fj,(C) and this causes R CUEp 7^ 
0. On the other hand if the perturbers increase i-i(A) without 
changing considerably the magnification of the other images, 
this will enhance the sum fi(A) + fi(C) pushing it closer to 
^i(-B), giving a smaller R C us P (0.07 in the case of figure 13). 

Figure 14 shows the effects of extra haloes on the R C ua P 



Figure 13. One perturbed configuration where iJcusp is less than 
the unperturbed value (0.091 vs 0.07). The solid circles show the 
position of the subhaloes, the point size is proportional to their 
mass. 



value for Config2. For all the three extra-haloes mass ranges 
considered, the modifications in the cusp relation are very 
small and the value of R cusp is not very sensitive to the total 
number of extra-haloes we generated (results for different 
number of subhaloes are shown by different curves, which 
are almost overlapping in the various panels). 

For this configuration we have performed one more test: 
we have modified by hand the fluxes of the three images in 
order to obtain an high value of i? C us P (> 0.37). Then we 
used the gravlens software to find positions and masses of 
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Figure 14. Probability distribution of R-variation considering 
Extra Haloes with different masses and distances from the pri- 
mary lens (Config2). Bottom right panel: comparison between 
the effects of extra haloes and subhaloes. 



Figure 15. Two subhaloes configurations with high R C uep value. 
The first one is indicated by solid squares and gives R C usp = 0.387. 
The second (one indicated by open circles) gives i? C usp = 0.274. 
In both cases M su i, Ri 10 6 . 



two subhaloes (with masses 5 x 10 5 < M < 5 x 1O 6 M ) 
with the constrain of simultaneously reproducing positions 
and fluxes of our modified images. We have found two con- 
figurations for subhaloes that are shown in figure 15. In the 
first case (solid squares) the mass of the substructures is 
roughly the same (w 2.0 x 10 6 ), they are close to the exter- 
nal images with a distance of 0.08 arcsec and the perturbed 
value of -Rcusp is 0.387. In the second case (open circles) the 
mass of the subhalo close to the central image is 1.6 x 1O 5 M0 
with a distance of 0.06 arcsec, while 7.8xlO 5 M is the mass 
of the second one that is far away from the central image, 
for a perturbed value of R C us P equal to 0.372. 

The aim of this test is to show that there is nearly 
always the possibility to explain an anomalous flux ratio us- 
ing subhaloes, but that their positions and masses must be 
tuned in a very precise way (i.e distances between images 
and subhaloes must be less than 0.08 arcsec). Most impor- 
tant, our Monte Carlo simulations show that the probability 
of obtaining such a fine tuning is very low. As a consequence 
we conclude that an explanation for the high number of ob- 
served anomalous flux ratios in lensed QSOs based on the 
presence of subhaloes in the mass range we have tested is 
very unlikely. 

4.4 Configuration 3 

In the last cusp configuration analyzed (Config3) the R-cusp 
relation is not completely satisfied even in the unperturbed 
case (-Rcusp = 0.243, fig 16). Figure 17 shows the R-cusp 
probability distribution for Config2 and Config3 normalized 
to the unperturbed value. There are no appreciable differ- 
ences between the two configurations. The probability dis- 
tributions have almost the same width and the same maxi- 
mum value (f« 0.35) and both are centered on the respective 
unperturbed value. This means that the ability of substruc- 
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Figure 16. Third unperturbed lens configuration: Config3 
(R=0.243) 



tures to modify the -Rcusp relation is nearly independent on 
the original value of the cusp ratio. 



5 FOLD RELATION 

For sake of completeness we have also considered a fold 
case similar to the configuration of the well known system 
PG1115 (Impey et al. 1998). In this case we define the fold 
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Figure 18. Unperturbed fold lens configuration (similar to 
PG1115, Rfoia = 0.92) 



relation as the ratio between the magnification of the closest 
pair on opposite sides of the critical curve (cfr fig 18): 



Rfo. 



fold 



= 1. 



(4) 



By applying our procedure to this test fold case we found 
that it is easier to modify the value of the unperturbed Rfoid 
with respect to the one of R C us P , as clearly shown in fig 19. In 
our computational procedure the position of the source can 
be changed by the gravlens code, in order to minimize the x 2 
of the lens configuration. While this does not affect the i? C usp 
relation it can be important in the fold case: as pointed out 
by Keeton, Gaudi and Petters (2005) the degree to which 
Rfoid can differ from one for realistic smooth lenses depends, 
in addition to the angular structure of the lens potential, not 
only on the distance of the source from the caustic but also 
on its location along the caustic itself. So the the values we 
got for Rfoid are due both to the effects of subhaloes and 
both to the shift along the caustics of the source. 

Thus, it is not possible to conclude from Rfoid alone 
whether observed flux ratios are anomalous or not. 



6 DISCUSSION & CONCLUSIONS 

Quasars that are being gravitationally lensed into multiple 
images have recently been used to place limits on the sur- 
face density of cold dark matter subhaloes (Mao & Schnei- 
der 1998, Metcalf & Madau 2001, Chiba 2002, Metcalf 
& Zhao 2002, Dalai & Kochanek 2002, Chen Kravtsov & 
Keeton 2003, Bradac et al. 2004, Mao et al. 2004). Small 
mass clumps that happen to lie near the images affect the 
observed magnification ratios. The question arises as to 
whether these observations are compatible with distortions 
expected to occur from dark matter substructures and satel- 
lite galaxies within the ACDM model. Recent results based 
on numerical N-body (Amara et al. 2004, Rozo et al. 2005) 
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Figure 19. Probability distribution of Rfoid f° r a different num- 
ber of substructures within an area of 3 arcsec 2 (5x 10 5 < M su t, < 
1O 7 M ). 



and hydro simulations (Maccio et al. 2005) have shown that 
it is hard to reconcile the observed high number of cups re- 
lation violation with the total amount of substructures pre- 
dicted by the ACDMmodel. These studies were limited by 
the present achievable numerical resolution that permits to 
resolve DM haloes down to masses w 1O 7 M . 

In this work we have quantified the effects of smaller 
mass clumps (10 5 — 10 7 Mq) on the observed violation of the 
R cusp relation. We employed results from N-body simula- 
tions to estimate the expected number of subhaloes in this 
low mass range. Due to the small mass of the perturbers 
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we have restricted our analysis only to those close (in 2D) 
to the images positions. For the mass range inspected in 
this work and for the typical distance between images (few 
arcsecs) this leads to a number of perturbers « 6. All the 
subhaloes are modeled as NFW spheres and we have gener- 
ated more than 10 5 different lensing configurations, varying 
masses, positions and number of subhaloes. 

The main finding of our work is that on a statistical ba- 
sis this class of perturbers is not able to modify consistently 
the unperturbed R-cusp relation. Values of i? C usp in the ob- 
served range (« 0.25) are obtained in only less than 1% of 
the analyzed systems. The ability of subhaloes in modifying 
the unperturbed value of R CUS p is found to be independent 
from the value of R C us P itself. 

These results are not in contradiction with the ones in 
the literature (Keeton et al. 2003, and more recently Mi- 
randa & Jetzer 2005 and references therein). As shown in 
figure f5 it is possible to use subhaloes in the mass range 
f 5 — f0 7 Mo to obtain high values of R C us P case by case, but 
a tight fine tuning between the location of the images and 
masses/positions of the perturbers is needed . 

In addition, we have also considered the impact of mas- 
sive haloes placed outside the primary lens (from groups of 
galaxies to a close cluster) by modeling them in the same 
way of the subhaloes. Our simulations show, as expected, 
that their contribution in modifying the R cusp relation is 
tiny and almost negligible with respect to the effect of sub- 
haloes. 

Results from this work together with results from nu- 
merical simulations seem to be in disagreement with the 
standard picture which explains the anomalous flux ratio 
by means of dark matter satellites. Interestingly, while on 
dwarf galaxy scale there is an excess of dark matter sub- 
haloes with respect to visible satellites, we have shown that 
the predicted level of substructures on smaller scales is not 
sufficient to explain the observed level of violation in the 
cusp relation. 

Possible solutions to this problem can reside in mi- 
crolensing for some of the lensing systems observed in the 
optical band (Metcalf 2005, Keeton et al 2005), or in the 
presence of haloes lying along the line of sight between the 
lens and the observer (Chen et al. 2003, Metcalf 2005), al- 
though total effect of this kind of perturbers is not yet clear. 
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